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Application is hereby made for continued step-funding support for the 
scientific programs in lunar research at the Lunar and Planetary Laboratory. 

A. Lunar Photographs and Maps Collection 

This has been augmented by the acquisition of Apollo 17 metric camera 
photography, also of the published sheets of the DMA Lunar Topographic Ortho- 
photomap series. Prospects for housing the entire collection of lunar and 
planetary imagery and maps in one secure location are now promising. It is 
hoped to obtain a computerized microfiche reader to simplify retrieval operations 
Very recently, an offer has been received from the new Main Library of the 
University, which is nearly completed and is situated close to LPL, to house 
this collection in a secure location there with librarian services provided 
gratis. However, the pros and cons of this offer will be examined very closely 
before any commitment is made. 

B. Personnel 

Lydia Chen and Sheryl Vrba have continued the crater cataloguing, 
measurement, and supporting card-punching operations under the immediate 
supervision of Leif Andersson, who has also been responsible for the thorough 
checking of all catalogued data, and the production of the important and 
fundamental statistics given later in this application. He has also re-written 
the crater-depth program (copy of method appended), prepared a program which 
produces coordinate grids for any Orbiter photograph, and written several 
other computer programs for processing the crater data. He obtained his PhD 
degree from Indiana University in September 1974. 

Bruce Cordell, a graduate student in the Department of Planetary Sciences, 
is doing theoretical work relating to the modification of the diameter/fre- 
quency distribution of craters by basin ejecta, and is comparing the result- 
ing expressions with crater counts from the new catalogue. 

Ralph Turner has completed all work on the model of Schroter's Valley, 

and has provided us with a fully illustrated report which only needs editing 
prior to publication. 
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C. Scientific Programs 

1 . Crater Measuring and Statistical Program 

General Status. It has been decided to publish the lunar crater 
catalogue in two volumes: one containing the earthside measurements and the 

other containing the farside measurements. These volumes will appear as Con- 
tributions of the LPL, since the material is too bulky for publication in the 
standard journal literature. 

To date measurements of over 17,000 cr aters have been compl eted _ . 

The entire earthside hemisphere is almost ready for publication and is expected 
to be in press by early summer. Only checking of a small part 
of the limb regions and the addition of a new piece of information (presence of 
a ray system) need to be completed. The number of craters measured on the 
earthside is about 11,800. The positions of about 1,000 craters very near the 
limb as listed in the old catalogue have been found to be seriously in error. 
These positions are currently being corrected. The outline maps and indices 
for the earthside are almost ready and will require only a little additional 

time to complete. 

The previously used depth calculation procedure has been replaced 
by a more exact one (Appendix I), prepared by Leif Andersson, and the photo- 
graphic supporting data used in these calculations have been updated. Depths 
for all craters with measured shadows are being recalculated using the new 
procedure and the updated photo data. Therefore, the catalogue of crater depths 
published by Arthur (1974), which was based on measurements made at LPL, does 

not give the best values available for these craters. 

About 40/S of the farside has been measured and most of these 
craters (over 5,000) are already on punched cards. Since much of the best 
available photography of the farside is at a rather small scale, only selected _ 
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parts (about 50%) will be measured to the same nominal limit of completeness 
(7 km) as the earthside; the rest will be complete to 10 km diameter. When 
completed, about 13,000 craters will have been measured on the farside. We 
had intended to use the farside coordinate grids from ACIC for crater positions. 
However, these have not been forthcoming and therefore Leif Andersson has 
prepared a computer program which generates these grids using the updated photo 
data. The program is very flexible and is also used to plot crater diameters 
on an overlay for any given photograph. This provides an extremely accurate 
and efficient method of checking crater positions and diameters. 

Preliminary Results . As stated in previous reports, crater depths 
based on the new measurements are significantly greater (often by 50-100%) 
than those calculated from earthbased measurements. Confirmation of these 
depths has recently been made by Pike (1974) using photogrammetric methods on 
limited numbers of craters contained on the Apollo metric camera photographs. 

As shown by Figures 1-5, the depth-diameter ratio for craters less than about 
25 km in diameter is fairly constant and well defined for each age class (Class 
1 is youngest, 5 oldest), indicating that the age classification is basically 
correct. 

The. diameter/frequency distribution of farside ^craters measured to 
date suggests that at leastthe central farside has a significantly higher 
crater density (~1.5) than typical southern highlands on the earthside (Figure 
6). This is evident for all crater sizes of Classes 1 through 3. This suggests 
that the paucity of craters in this area of the earthside highland is not due to 
nearby basin formation since this would have preferentially obliterated the 
smaller craters with respect to the larger. The investigation of the number 
density of craters in selected areas of the highlands will be very valuable 
in evaluating possible causes of these differences. 
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Fragmentation theory predicts that the log frequency versus log diameter 
plot for impact craters should show a -2 slope. The most widely referenced 
earlier crater counts are those of Hartmann (1966, 1967) which seemed to 
confirm the -2 slope. These counts in highland areas were based on the old 
LPL crater catalogue and on early Lunar Orbiter photos of the farside. Further- 
more, the counts were made in selected areas of the highlands (termed "Pure. 
Uplands") which show high crater densities and were considered to be equili- 
brium surfaces. The counts on the frontside show a slight downward bend 
below about 60 km, but when combined with the farside data the best fit seemed 
to be a straight line with a -2 slope. Hartmann then compared this -2 slope 
with crater counts for the heavily cratered provinces on Mars and reached rather 
fundamental conclusions concerning an episode of atmospheric erosion on Mars. 

Figures 7 and 8 are diameter /frequency plots of the Hartmann counts for 
the Moon and Mars, Gault's counts for the cratered terrain on Mercury, and the 

new LPL counts for the lunar highlands based on the entire earthside and 40% 

% „ 

of the farside. Our lunar data are much more complete than Hartmann's and 
represent a large homogeneous sample of highlands rather than selected areas. 
Figure 7 is the traditional diameter/frequency plot using logarithmic incre- 
ments, while Figure 8 is a log-log plot of diameter versus relative crater 
density in craters per diameter squared of the different size increments. 

In Figure 8 these curves are referenced to Hartmann's -2 slope which is shown 
as a horizontal line. This type of plot shows much more graphically and realis- 
tically the differences and similarities between the curves, and significant 
departures from a -2 slope. The new lunar highland curve shows several new 
and far-reaching results. 

Within the diameter limits of the data (7-450 km) the lunar highland 
counts do not conform to a -2 slope. The peak density occurs at about 80 km 
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and shows a significant decrease in density at diameters both below and above 
this value. If Hartmann's data points for the nearside and farside are 
examined separately, this same trend is seen in both sets of data. As pre- 
viously mentioned, our data suggest that the crater density in almost all size 
ranges is greater on the farside highlands than the earthside highlands. 
Furthermore, our farside highlands data alone (Figure 6) clearly show a 
decrease in relative density below about 70 km in diameter when referenced 
to a -2 slope. We believe that when Hartmann combined his nearside data with 
the higher density farside, the general crater density was raised and the 
downturn of the curve was masked to give a "best fit" of all data to a -2 
slope. The decrease in. density below about 80 km in diameter shown on our 
curve may be due to obliteration of craters caused by basin formation and/or 
an early period of highland volcanism, or it may represent the size distri- 
bution of the impacting bodies. The apparent paucity of craters in the size 
range about 80-450 km may be due to one or more of the following: (1) exten- 

sive slumping in the larger craters and multinring basins thereby causing a 
larger diameter than the true excavation diameter, (2) a different population 
of large impacting bodies, or (3) the influx of large bodies being confined 
to an early history and then subsequently obliterated by later impacts and 
volcanism. When other statistical parameters in the catalogue are taken into 
account it may be possible to eliminate several of these alternatives. 

Very significant is the fact that all three bodies (Moon, Mars and 
Mercury) show the same general shape of the crater diameter/freque ncy distri- 
buti on . This indicates that the general cratering and subsequent modification 
histories are similar. Furthermore, the similarity between the lunar and 
martian curves indicates that the decrease in crater density below 60 km on 
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Mars is not solely the result of an episode of atmospheric erosion as 
proposed bv Hartmann. As seen in Figure 8 the crater density of Mars is 
at least a factor of 2 less than that of the Moon at all sizes. If this 
crater density difference is taken into account, then the loss of craters 
between about 50 and 7 km diameter due to atmospheric erosion on Mars is 
only a factor of 1-2. This is considerably less than the factor of 2-18 for 
the same diameter range given by Hartmann's lunar and martian curves. These 
preliminary results are being prepared for publication. 

It is clear from the preceding discussion that crater parameters for 
the entire Moon are required to adequately portray its cratering and sub- 
sequent modification history. Furthermore, these data are essential before 
reliable comparisons can be made with other planetary bodies. 

II 

2. Schroter's Valley Model 

The contouring of the entire model by computer techniques was 
completed late in 1974, and Figure 9 is a much reduced copy of the assembled 
mosaic, presented in a Mercator projection. Absolute altitudes correspond 
to the best fit with data from DMA chart LT038B3, which overlaps a small 
area of the model. The depth of the Valley in this area is approximately 
600 m both from the model and the DMA chart; this gives one some confidence 
in the general modeling technique used. 
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Cross sections in N-S and E-W directions, as well as height profiles 
along the Valley are given. These show very graphically the gradually de- 
creasing downwards slope from the "Cobra Head" to the "tail", and the general 
slopes at specific locations (Figure 10). Arguments have been advanced for 
supposing that the narrow ends of sinuous rilles may be higher, thus resembling 
terrestrial rivers, with the flowing material draining down the crater at the 
broad ends. The model shows that this is certainly incorrect for Schroter's 
Valley, and by analogy, for other sinuous rilles. 

3. 61 -inch Color Filter Photography . 

Mr. Whitaker was able to pursue this project actively for a period 
late in 1974, allowing relative calibration of the original films obtained with 
the 61-inch telescope. Methods of exposure and processing (with the LPL Transflo 
processor) of the intermediate positive transparencies were standardized to 
produce gamma (i.e. contrast) curves of the requisite shape and slope for any of 
the film pairs. It soon became apparent that with absolute calibration, the o 
original films could yield reliable color and albedo values in addition to the 
color boundaries . A preliminary attempt was made in late March to make such cali- 
brations, using a Pritchard photometer (kindly made available by Dr. L. Dunkelman 
of GSFC) attached to the 61 -inch reflector. The results look very promising, 
and may already be suitable to calibrate the original images. 

With both albedo and relative color data available, it will be possible 
to prepare a two dimensional plot which will segregate different units much more 
effectively than by color alone. 

Figure 11 is a color-difference photograph of the Aristarchus region, 
showing the degree of resolution obtainable by combining photos. A comparison 
between this image and several other data sets (radar, IR, geology, etc.) was made 
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late in 1974 in conjunction with members of the S217 and S222 Lunar Synthesis 
Teams. The results were presented at the Sixth Lunar Science Conference, by 
S. Zisk and will be published in due course. 

This project was funded separately last year (NSG 7014), but expected 
expenditures to complete project are sufficiently small to include them in this 
proposal. 



REFERENCES 


Hartmann, W. K., Early lunar cratering, Icarus , vol. 5, 406-418, 1966. 

Hartmann, W. K., Martian cratering, Icarus , vol .5,565-576, 1966. 

Hartmann, W. K. , Lunar crater counts. II: Three lunar surface type- 
areas, Comm.LPL , vol. 6, no. 81, 39-41, 1967. 

Hartmann, W. K. , Martian cratering, 4, Mariner 9 initial analysis 
of cratering chronology, J . Geophys . Res . , vol. 78, no. 20, 4096- 
4116, 1973. 

Murray, B. , et al.. Mercury's surface: Preliminary description and 
interpretation from Mariner 10 pictures. Science , 185, 12 July 
1974. 



Cenfrrxtf E<*rtk$ide 

II 27 craters 







5 


10 


20 


50 


100 DIAMETER 




C*£ntr*s/ Earth £ide 

^9 



A-.oyD 


Fi 3 . 3 

I-+ — 


ion !']tamf:tf:r 












I 


10 


100 


1000 




RELATIVE DENSITY (Craters/ I0 3 D 2 Km 2 ) 


Pure Lunar Highlands 
(Hartmann) 

+ Earth based 
□ Orbiter farside 


Lunar Highlands 
(new LPL data) 


Mercury — 
"Highlands" 
(Gault) 


Mars: Heavily Cratered Provinces 

(Hartmann) 


Figure 8 


DIAMETER (KM) 


IOO 




cgroyo 



s*ktw f kYv. 




^J"20'rir- <-, SJVOw 


&ylf 


WMIU 






war 




fMr 

n 




sv 






,i AUM 


mA « cs \ 


sroo’tv 


Sl°20w 


SfOOw 

center oFchord plane i - 2$"00'n 


SO'+Ow 50"2ffw\ 




♦r qr 




r; 


tn<-i\\X-X‘ TXVV I 


SCHROTER'S VALLEY 


Ca)W 
Pf ~Mw 

h /.^\a 


S OT<?w 


K»n r> .5 ><> 


KXW 5000 


7 y; 



^ f Vt n) fvrv 
£ 1 ( \ y M y V 




.; 1 VYS 

5**0 W: sA\>’w 


i 


Measurements from Relief maid— RfJurner iyjf 

Ho. 9 


LO^' 

4 Yv 

SO'Oflw 



Krtt 

i o: $ 




i rnh 

IWQw 


A9"20'w 


lumvmm 


HWw 


^^OOt/y 


$TtA 1/V4 










Color-difference photograph of Aristarchus-Wood's Spot region of Moon 
(red minus UV). Redder areas appear lighter. 





Depth calculations for lunar craters from Orbiter photos. 

(current method by L. Andersson) 

The input data are: 

for the photograph: spacecraft altitude a 

'* long i tude 
" latitude 

11 ti 1 1 angle T 

" tilt azimuth 

Sun's colongitude 
11 latitude B 

o 

for the crater: coordinates of shadow casting point >1 

shadow length on photograph I' £see below for unit used^ 


Cartesian coordinates are in units of the Moon's mean radius R - 1738.1 km 
Origin is at the Moon's center. £ axis poi nts. towards East, n axis 
towards North, and £ axis towards Earth. Shadow casting points are 
assumed to be on Moon's mean surface so that £ 2 + n 2 + S 2 =1* 
distances are in units of R except where specifically noted. The 
following subscripts will be used: 


subscript 

o 

s 

n 

P 

a 

none 


referring to 

Sun or subsolar point 

spacecraft 

" nadir 

principal point 
camera axis 

crater (shadow casting point) 


In the figures, a bar is used to indicate a complement angle (e.g. 

= 90° - B). The symbol N denotes the North pole, n the spacecraft 
nadir, s the spacecraft, 0 the center of the Moon, P the principal 
ground point, o the subsolar point, and C the crater (shadow casting 
point) . 
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The following quantities are obtained directly from the input data 

p = 1 + a/R 

£ = p sin X cos B 

n n n 

H = p s i n B 
n n 

£ = p cos X cos B„ 

n n n 


A - 90° 

- X 

o 

o 

£ = sin 

X cos B 

o 

O ' 

n - sin 

B 

o 

O 

C = cos 

X cos B 

o 

o 


C - 1 - " n 2 1 


X “ arctg £/C 
B = a res in n 

Loc ation of Sun with respect to crater (Fig. l).L 

Sun’s elevation angle = E q 
S un’s azimuth ~ A q 

sin E q - sin B sin B Q + cos B cos B Q cos (X-X q ) (cosine theorem on A CNo) 

-« 0 + nn o + CC 0 

sin A - cos B sin (X-X )/cos E n (sine theorem on A CNo) 

o o o o 

cos A - (sin B - Sin 3 sin E )/cos B cos E (cosine theorem on A CN 0 ) 

o o ° v 
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Location of spacecraft with respect to crater (Figs. 2a and 2b) : 

Selenocentr Ic angle between crater and spacecraft = g 
Spacecraft azimuth = An 
Spacecraft elevation angle ~ E^ 

distance from crater to spacecraft ® x [units of R] 


cos q = sin 8 sin B + cos B cos 6 cos U -A) (cosine theorem on A CNn) 

n n n 

■ • « « +‘ nn + CC 

^ n n n 

sin An = cos 6 sin (X -X)/sin g . (sine theorem on A CNn) 

n n 

n 6 cos g)/cos P sin g (cosine theorem on A CNn) 

(plane cosine theorem on A COs) 
(plane sine theorem on A COs) 


cos An = (sin (S - si 
n 


V - ■ » 

1 + p 2 - 2p cos g 
cos = p sin g/x 


Relations between crater, spacecraft, and camera axis (Figs. 3a and 3b); 


Let Q. be the point whose coordinates are the direction cosines of the 
camera axis (the positive direction of which is taken to fee from P to s) . 
Selenocentr ic angle between spacecraft and principal point = k 

Slant range (distance between spacecraft and p.p.) = S [units of R] 

angle between camera axis and position vector of crater = h 


azimuth of Q, at crater = A^ 


sin = ccs T sin B - sin T cos 0 cos A 7 (cosine theorem on A QNn) 
Q n n Z 

sin(X ~X_) = sin T sin A 7 /cos (sine theorem on A QNn) 

n Q. Z Q. 

cos(X -X_) = (cos T - sin 6 sin g n )/c.os 6 cos ' 

n Q n U n u 

(cosine theorem on A Q.Nn) 


X 

n 


( W 
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*Q = sin X Q cos B Q 


\ = sin B Q 

C Q = cos cos B q 


S = p cos T " p 2 sin 2 T 

sin k = S sin T 

sin = cos k sin 3 + sin k cos 3 cos A 7 

P n n/. 

sin (X„-X ) = sin k sin A 7 /cos B D 
P n Z P 


(plane cosine theorem on A POs) 

(cosine theorem on A PNn) 
(sine theorem on A PNn) 


cos(Xp _ X n ) = (cos k _ sin sm Pp/z^^ p n >- i p 

X p = A n + .<V X n ) 

£ p = sin X p cos B p 

n p = sin B p 

C p = cos X p cos B p 

cos h = sin B sin Bq + cos 3 cos cos (X^-X) 

+ n \ + CC Q 

sin = cos Bq sin(X^-X)/sin h 

cos Aq = (sin Bq " sin B cos h)/cos B sin h 


B sin Br,)/cos B cos B p (cosine theorem on A PNn) 


(cosine theorem on A CNQ.) 

(sine theorenTon A CNQ) 
(cosine theorem on A CNQ) 
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Linear distance between P and C — Z C un >ts of R] 

Angle at spacecraft between P and C = a (Fig- *0 

Angle between sun direction and camera axis = x (Fig. 5) • 

i z = U -£ p ) 2 + (n-n p ) 2 + (?-Cp ) 2 

cos a - '(x 2 + S 2 - Si 2 )/ 2xS (plane cosine theorem on A CPs) 

cos x = sin & sin + cosg cosg n cos(X -A ) (cosine theorem on A QNo) 
o y. o y. *1 v 

Now consider Fig, 6. The Z^axis is the camera axis. The XY-plane, 
which we will call the “object plane", contains the point C (the 
shadow-casting point). The end (bottom) of the shadow is at K; the 
length of the shadow hypotenuse is 0 . The projection of the shadow 
on the object plane is CL. The angle in the plane, at the crater, 
between this projection and JC (the projection of the spacecraft-crater 
line) is denoted by H. Finally, M is the intersection of JL and sK, 

The line CM, of length Z, is the “apparent shadow length", which is 
in the same ratio to the distance spacecraft-object plane as is the 
length of the shadow image in the film plane to the focal length of 
the camera. 

Let the distance LJ be = t, LM = p, and the angle JLC = y. Then: 

t 2 = a 2 sin 2 x + x 2 sin 2 a + 2o X sin x sin a cos H (plane cosine theorem 

on A JCL) 

sin v = X si n H s ' n a/t (plane sine theorem on A JCL) 

p _ a cos x (A JsM <s» A LMK) 

t-p x cos a 

E 2 = p 2 + o 2 sin 2 x - 2p0 sin x cos y (plane cosine theorem on A CLM) 





Combining these four equations, eliminating p, t, and y, yields after 
some manipulation: 


2 _ y 2 g 2 (cos 2 a sin 2 x + sin 2 a cos 2 x - 2 sin x cos x sin a cos a cos H) 
S (a cos x + x cos aP ' " 


^ a m E X cos a where G = Vsin 2 (x+a) - • +C ° S - - sin' 2x sin 2cJ 

XG ~ E cos x v 2 


To obtain cr in ki 1 ometers , x and E must be converted to this unit. 


X km = 1738.1X 

(Since in Fig. 6 Js may be 
thought of as the camera's focal 
length, f, and £ as the shadow 
length in the film plane, S') 

Here a is the ratio of the scale of the spacecraft film to the scale 
of the prints used for the shadow measurements, a — 0.1380 is the 
adopted value. 


Thus 

• „ 1738. lx 

km G f 

cos a aE 1 


Z km = *738.1 x a cos « 


1 



~ 1 ~ 


Let the altitude of the shadow casting point over the mean surface 
of the Moon be 6 (Fig. 7), the depth of the crater d, and the pre- 
viously derived elevation angle of the Sun E q . Then 


(R + 6- d) 2 = (R + S) 2 + cr 2 - 2 (R +S)c? sin E q 


d = (R +"S ) 



o sin E 
o 

R + 6 



) 


(plane cosine theorem) 


Expanding in a series one obtains 

d = a sin E - cos 2 E + cos 2 

o 2R o ^2 


E (o sin E + 26) + . . . . 
o o 


Since <5 is generally unknown, but small (of the order of d) , it is 
appropriate to retain two and only two terms of the expansion: 

, • r- o ' 2 o r 

d = a sin E - -ztt cos* E 
o 2R o 





Ftq, 1 

cj 



n 3 . 1 a 
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